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Powerful computational fluid dynamics (CFD) tools have emerged that appear to offer 
significant benefits as an adjunct to the experimental methods used by the stability and 
control community to predict aerodynamic parameters. The decreasing costs for and 
increasing availability of computing hours are making these applications increasingly 
viable as time goes on and the cost of computing continues to drop. This paper summarizes 
the efforts of four organizations to utilize high-end computational fluid dynamics (CFD) 
tools to address the challenges of the stability and control arena. General motivation and 
the backdrop for these efforts will be summarized as well as examples of current 
applications. 


Nomenclature 

C D = drag coefficient 

Cl = lift coefficient 

C La = derivative of lift coefficient with a 

Ci = rolling moment coefficient 

C m = pitching moment coefficient 

C n = yawing moment coefficient 

C n/3 = directional stability 

C p = pressure coefficient 

C Y = side force coefficient 

c = mean aerodynamic chord, ft. 

k = nondimensional frequency, coc/2U 
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M x = free-stream Mach number 

Re c = Reynolds number based on mean aerodynamic chord 

£/oo = free-stream velocity, ft/sec 

a = angle of attack, deg 

f = angle of sideslip, deg 

8 S = horizontal tail deflection, deg 

co = dimensional frequency, 1/sec 


I. Introduction 

R egardless of expectations for range, cruise Mach number, or payload, if an airplane does not have acceptable 
stability and control (S&C) attributes, it will not go into production. Although the state of the art in prediction 
of aerodynamic S&C parameters has vastly improved and the great majority of S&C deficiencies identified and 
corrected during developmental wind tunnel testing, unsatisfactory S&C characteristics are occasionally not 
identified before going to flight test. When these shortcomings do surface during flight test, however, they have 
often resulted in costly and inefficient problem resolution efforts with varying degrees of success. Many of the 
problems have been directly related to extremely complex aerodynamic phenomena and the difficulty in 
predicting and resolving such issues. The critical challenge for aerodynamicists is to reliably predict, and if 
necessary improve, all critical S&C characteristics before getting to the flight test stage of development. 

The wind tunnel has traditionally served as the primary tool to provide aerodynamic inputs for simulation 
data bases and to predict S&C characteristics. Unfortunately, shortcomings exist in wind tunnel testing methods 
that can result in serious errors in the predicted S&C characteristics. For example, most wind tunnel testing for 
S&C is conducted at subscale values of Reynolds number. Consequently, large effects of Reynolds number can 
occur with certain configurations because of the dependence of boundary-layer separation on Reynolds number. 
In the case of dynamic derivative testing, both low Reynolds numbers and low Mach numbers are typical due to 
facility constraints. As summarized by Wolowicz, Bowman, and Gilbert 1 , testing at reduced values of Reynolds 
and Mach numbers poses significant challenges. Model fidelity and model support interference can also be 
significant issues for wind tunnel testing. 

It is the belief of the authors that the next significant improvement in the state of the art for predicting the 
S&C characteristics of a new vehicle is underway at this time with the application of powerful computational fluid 
dynamics (CFD) tools. While the current time- accurate, Navier-Stokes solutions to typical S&C problems are still 
too resource intensive for widespread use, these costs are expected to continue to fall, at least for the next decade, 
according to Moore's Law. 2 Assessment and validation efforts are now timely so that validated tools will be 
available when computational costs have dropped to the point where such solutions would be reasonably 
affordable in a typical S&C setting. 

Applying high-end CFD codes with technology greater or equal to Reynolds-averaged Navier Stokes (RANS) 
to specific areas of S&C interest before first flight can help focus the wind tunnel program and provide improved 
understanding of the flow physics encountered. Developing, assessing, and maturing appropriate CFD codes for 
use in S&C analyses is the goal of the Computational Methods for Stability and Control (COMSAC) project at 
NASA Langley and a goal of the various organizations that coauthor this paper. The purpose of the current report 
is to present the reasons why this application of high-end computational tools is timely, to highlight the efforts 
that are underway in each organization, and to show how these efforts can synergistically accelerate the progress 
of validating the high-end codes for both simple and complex stability and control challenges. 


II. Historical Examples of Shortcomings in Predictions for Stability and Control 

Historically, many civil and military aircraft development programs have encountered critical S&C 
deficiencies during early stages of flight test. Despite thousands of hours of wind tunnel testing, these surprises 
have occurred across the speed range from takeoff/landing to cruise, with especially challenging encounters at the 
fringes of the flight envelope, where separated flows dominate. 

Examples of a number of high performance military aircraft that encountered transonic, uncommanded lateral 
motions during flight test are summarized in figure 1 . 3 The common experience of the aircraft in figure 1 is that 
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critical S&C deficiencies were uncovered in flight test. Consequently, fixes had to be quickly found under less 
than ideal conditions in a "cut and try” mode without the general benefit of flow diagnostics. In many cases, 
adding vortex generators or fences rectified the problem while, in other cases, performance limitations were 
simply accepted. Reliably predicting such unacceptable S&C characteristics before building flight hardware 
could result in more efficient configuration optimization and could avoid costly program delays. 

In addition to the foregoing high-performance military aircraft, a broader spectrum of aircraft that 
encountered S&C surprises in flight test is shown in figure 2. The examples shown include the Boeing 767-400, 
where the addition of raked wing tips altered the stall break and necessitated the development of a new vortilon 
pattern. 3 4 The Lockheed Martin C-130J experienced a wing drop, or roll-off, due to propeller-induced effects, 
which impacted product development of that aircraft. The experience with the Boeing 111 involved missed 
predictions concerning aft loading on the main wing and effectiveness of the horizontal tail in flight. 4 The Lear 23 
development program had the challenge of fixing an abrupt wing stall associated with leading-edge separation on 
approach. The NASA experience with the loss of the first X-43A stack configuration appeared to be the result of 
a cumulative sum of uncertainties in many of the predicted stability and control characteristics. Finally, the 
Boeing 737 Next Generation aircraft demonstrated a new, and unexpected, sensitivity of rigging of the inboard 
trailing-edge flaps that was not apparent in the older models. 4 These revelations from flight test manifested 
themselves as increased time and cost of development testing, acceptance of undesirable performance limitations, 
or limitations in future growth potential. 

On the basis of this past experience, it is clear that vehicles across the spectrum from high-performance 
military to commercial transports to business jets have all suffered from S&C deficiencies first realized in flight 
test. Even more sobering is the fact that nearly all of the vehicles in figure 2 could be considered derivatives of 
very similar earlier versions of the configuration. As the industry moves away from derivative aircraft to more 
revolutionary vehicles, whether blended wing bodies, uninhabited vehicles, or vehicles with active control 
effectors, either more accurate predictive tools must be developed or even more ’’surprises” will undoubtedly be 
encountered. 


III. Bridging the Cultural Divide 

It is now hoped that CFD tools can be brought to bear on S&C issues and complement the traditional wind 
tunnel data sources. However, for this next step to become a practical reality, improved communications between 
the S&C and computational communities must occur, even though they have very different backgrounds. Some 
of the cultural differences between these two communities are summarized in a tongue-in-cheek comparison in 
figure 3. During the past 10 to 15 years, the CFD community has generally been focused on performance 
predictions associated with attached flows. It has been able to focus on flow physics, which is a beneficial by- 
product of the computations. It’s primary concern about flow separation has involved predicting the value of a at 
which it might begin. The focus of the performance computations was typically forces such as lift and drag for 
one or two design conditions involving symmetric and steady flight with attached flow. Accuracies were 
expected to be good, and analysts sought to optimize the configuration for performance. Knowledge of S&C was 
not required. 

In contrast, the S&C community has had to worry about integrated aerodynamic forces and moments in all 
six degrees of freedom including the effects of aircraft motion. Since the S&C community is responsible for the 
vehicle over its entire operating envelope, it has had to be concerned about both attached flow conditions as well 
as massively separated flow conditions at values of a up to and beyond stall and for large values of (3. Symmetric 
flight is important but is only a portion of the flight envelope. Whereas high accuracy is appropriate and needed 
for performance predictions, the S&C community would sometimes be reassured just to know the algebraic sign 
of a parameter— for example, to know if the configuration is dynamically damped or propelling in pitch, roll, and 
yaw. To improve a configuration or to resolve unacceptable deficiencies, there have been very few guidelines and 
only minimal flow diagnostics available. Consequently, cut and try has been the typical approach, which 
depended heavily on the practitioner’s experience base with similar configurations. Little or no knowledge of 
CFD has been required for the majority of these engineers. 

If CFD tools are to be applied to S&C problems, both the S&C and the CFD communities will have to 
actively be involved. It will be critical that the S&C community define levels of accuracy needed for both static 
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and dynamic derivatives for a given vehicle. If accuracies can be relaxed relative to performance computations, 
then computational resources and, therefore, cost can be minimized assuming that the CFD community can 
quantify its uncertainty levels. Also, the S&C community must help focus the application of the computational 
efforts. While wind-tunnel testing is designed to provide large data bases that can cover most of the flight 
envelope, high-end computations are resource intensive and can currently only be strategically applied. Finally, 
technical trust must be developed between the two communities. The CFD community will have to demonstrate 
that it can predict results in a timely way and without prior knowledge of the answer if its results and tools are to 
be accepted by the S&C community for any role beyond providing flow diagnostics. Also, the S&C community 
will need to become educated and conversant with CFD tools designed to be user friendly. 


IV. Recognizing the Need for CFD Tools 

During November of 2002, a NASA/DoD/Industry Flight Prediction Workshop sponsored by NASA and 
DoD was held in Williamsburg, Virginia. The invitation-only meeting was attended by 85 experts who shared 
critical experiences and involving the state-of-the-art in aerodynamic flight prediction. During break-out sessions, 
stability and control deficiencies were highlighted as high priority and the lack of robust, accurate computational 
tools as an adjunct to wind tunnel testing were also cited. These conclusions were used to help launch the current 
COMSAC effort. 

In September of 2003, a follow-up Computational Methods for Stability and Control (COMSAC) Symposium 
sponsored by NASA was held in Hampton, Virginia with approximately 100 attendees. Thirty- four presentations 
were given including keynote and overview presentations from government and industry representatives as well as 
technical papers outlining the state of the art from government, industry, and academia. 5,6 Both computational and 
S&C perspectives were presented. The challenge of working cooperatively while protecting proprietary concerns 
was discussed. 

During the course of the COMSAC Symposium, it became clear that a high level of interest exists in pursuing 
the application of high-end codes to S&C challenges. Coupled with the already existing shortfalls in wind tunnel 
capabilities is the growing realization that wind tunnel availability may be more limited in the future as a number 
of existing facilities are being considered for closure. All these factors were stimulating interest in accelerating 
the application of CFD tools to S&C problems. It was also clear that many organizations are already actively 
pursuing interactions between these two disciplines. 

However, several concerns were also expressed by attendees and speakers at the COMSAC Symposium. The 
first concern was that CFD had been ’’over marketed” on many occasions in the past and that, understandably, the 
S&C community was skeptical of current capabilities in view of the extremely challenging separated flows 
encountered in its efforts. Also, the S&C community is aware that many current diagnostic computations are 
made with the code practitioner knowing in advance the answer and, therefore, able to optimize the code results 
through grid modifications or turbulence model selection. They wanted to see more ’’blind” predictions. The 
second major concern was that the high-end codes using RANS appear to produce incorrect answers for a large 
fraction of the S&C cases undertaken. 

In summary, the symposium brought together two different engineering cultures for in-depth discussion of 
opportunities, challenges, and experiences with CFD methods for S&C. Agreement existed in both cultures that 
applications of high-end computational tools to current S&C challenges for calibration and validation would be a 
milestone goal. Currently, it is obvious that this validation is very much in its infancy and that a considerable 
amount of work comparing computational predictions for many different configurations to actual flight, or wind 
tunnel, data will be required before these tools will be fully accepted by the S&C community. 

V. COMSAC Project Objective 

Based on the COMSAC Symposium and other NASA/industry/government interactions, the objective for the 
COMSAC effort might best be stated as ’’accelerating the application, validation, and focused development of 
CFD methodology to S&C aerodynamic predictions and analyses by coordinating industry/government efforts.” 
The benefits of this effort are envisioned to include (1) better understanding and control of critical flow physics, 
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(2) reduced and focused wind-tunnel and flight tests, and (3) risk reduction for future programs because of fewer 
surprises in flight test and certification. An informal COMSAC community was formed by NASA Langley, 
Naval Air Systems Command (NAVAIR), Boeing Commercial Airplanes, and the Air Force Research Laboratory 
(AFRL). While it is fully understood that many other industry and other government entities are involved in the 
pursuit of applying high-end computational codes to S&C problems, the above listed organizations are attempting, 
through informal information sharing, to accelerate the process. 

VI. Status of Efforts 

A. NASA Langley 

The Langley approach has been to apply several mature, state-of-the-art computational codes to the respective 
problems of interest to explore differences resulting from the respective implementations of algorithm, turbulence 
modeling, and other factors. The unconventional Blended Wing Body (BWB) configuration has been the focus of 
the effort because of previous NASA/Boeing research and current interest in this advanced configuration. While 
the BWB has undergone several configuration evolutions over the past decade, the configuration being analyzed 
for the present work is the 450- 1L configuration. 7,8 As seen in figure 4, this flying-wing type configuration has a 
relatively thick center section with conventional outer wing panels. The three-engine installation incorporates 
pylon-mounted nacelles although research has also been accomplished in the past for boundary- layer ingestion 
inlets as well 9 . 

The COMSAC BWB computational efforts have been underway since February of 2004 and are compared to 
Langley wind tunnel data which were intended to populate a simulation data base. Unfortunately, the existing 
experimental data sets, (particularly for the pitching moment coefficient) are all influenced in varying degrees by 
the presence of the model-mounting system. The photographs shown in figure 4 illustrate a .03 -scale model of the 
BWB mounted on two examples of the post supports used during the Langley 14- x 22-Foot Subsonic Tunnel 
(ST) and a 0.02-scale model of the BWB on a blade support in the Langley National Transonic Facility (NTF). 
The post mounts in the ST varied from bottom posts entering the model vertically through its belly region with 
different post diameters to a post entering the top, back of the model. All of the supports generated aerodynamic 
interference increments in the experimental data to some degree, as may be seen in Figure 5, which shows the 
variation in pitching moment measured in the NTF and on large and small posts in the 14- x 22-ft ST. Although 
the NTF data and the large-post 14- x 22-ft ST data are for the same nominal Reynolds number (3.7 million), the 
small-post 14- x 22-ft ST data are at a slightly lower Reynolds number (2.5 million). To help gauge the relative 
effects of Reynolds number vs. post-size, Figure 5 also shows the large-post 14- x 22-ft ST at a Reynolds number 
of 2.7 million. While there is a noticeable effect due to Reynolds number differences, which were created by 
changing tunnel velocity, the post size has a far greater effect at virtually all values of a. 

The knowledge that mounting support interference was an issue was identified during the early 14- x 22-Foot 
ST testing program and led to computational corrections for the large post experimental data by Boeing 
Huntington Beach. Huntington Beach also designed the blade support for the NTF test to minimize interference 
effects. The COMSAC effort at Langley is continuing the effort of Huntington Beach to quantify the impact of 
the various support systems used so that the wind tunnel data can more easily be used to validate the CFD codes. 

7. Static computations 

Before discussing the correlation of CFD and experimental results for the BWB configuration, a brief 
discussion of pertinent S&C interests might serve to provide the reader with important background. Shown in 
Figure 6 is one of the experimental pitching moment curves from figure 5 and can be used to illustrate three key 
regions for S&C. For the lower angles of attack near 0°, the pitching moments are generally linear with increasing 
AOA and the slope is negative, indicative of satisfactory longitudinal stability; such trends are usually exhibited 
for attached-flow conditions. Within the second AOA regime wing stall has resulted in shifting of the 
aerodynamic center of the configuration. This region is initiated by an abrupt longitudinal instability (pitchup) 
indicated by a reversal of slope to a positive trend and extends over the region of complex vortical and other 
separated flows. Finally, the third AOA regime occurs when the complete configuration stalls and the pitching 
moments return to a stable slope. The S&C analyst is concerned with the stability trends shown in each of the 
regimes — especially to the abruptness of instabilities and the proximity to the intended flight envelope as well as 
recovery from post-stall conditions. 
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For all of the Langley computations to be presented, no attempt to model the tunnel walls has been made and 
the results are considered to be appropriate for ’’free air." Furthermore, initial computations did not account for 
the presence of the mounting supports while later computations, as will be shown, did. Most attention has been 
focused on capturing the longitudinal pitching moment coefficient, as this is critical to characterize longitudinal 
stability and has historically been difficult to accurately predict with CFD. In addition to comparison with 
experimental data, most cases have been run with at least two CFD codes for code-to-code comparisons. While by 
no means a true validation process, having two different codes giving nominally the same result leads to 
additional confidence in the computations that might otherwise be absent; on the other hand, large differences in 
computed results would suggest the need for further investigation. 

With the understanding that the experimental data contain support interference while the initial CFD 
computations do not, the computational pitching moment results from the unstructured-grid solvers FUN3D 10 and 
COBALT 11 as well as the structured-grid solver PAB3D 12 are compared to experimental data in Figure 7. This 
figure shows a comparison of experimental data for Re c , Reynolds number based on mean aerodynamic chord, 
equal to 3.7 million, and for = 0.2. The FUN3D and COBALT results were generated using the Spalart- 
Allmaras 13 (SA) turbulence model, while the PAB3D results were generated using the Girimaji 14 explicit algebraic 
stress model. Additional static computations obtained with PAB3D and FUN3D, with comparison to experiment, 
can be found in reference 15. 

The results from COBALT tend to agree closer with the NTF data in terms of stability levels at lower values 
of a than the FUN3D predictions, but both capture the overall trends very well to the higher values of a. The 
results shown correspond to the finest meshes used in the two solvers; COBALT is cell-based, so the number of 
cells is the proper measure of the grid size, while FUN3D is node based, thus the number of nodes indicates the 
problem size. More variation is seen between the two unstructured solvers and the structured solver, which shows 
a much lower level of stability at lower values of a. It should be noted that the PAB3D results were for a coarser 
grid (700,000 cells) and with a different turbulence model. The data set from the 14- x 22-ft ST was taken with a 
large post support mounted through the belly of the model, and although displaced with a positive increment of 
moment, the trends (stability level at lower values of a and break points) agree with the NTF data. Since the NTF 
data utilized a lower surface blade support, which was designed to minimize support interference, it is 
encouraging that the CFD results fall closer to the NTF data than they do to the 14- x 22-ft ST. However, the 
large differences in the experimental values, due to different levels of support interference, preclude any real 
validation of the computational data without modeling the respective support systems. 

Because the experimentally measured dynamic data discussed in the next section were of necessity taken at 
very low dynamic pressures, a similar comparison of CFD and experimental static moment coefficients was 
repeated for Re c = 700,000, corresponding to the flow conditions for the dynamic tests. The grids used were the 
same as that for the higher Reynolds number, with the exception that the COBALT results shown were using a 
coarser mesh that was subsequently used for the majority of the dynamic pitch computations run with COBALT. 
This grid was chosen based on running three different grid resolutions at the higher Reynolds number (static) and 
by running two dynamic cases with the three grids. The experimental measurements were taken using a mounting 
system that allowed the models to be mounted either from the top or the bottom; both experimental data sets are 
shown in Figure 8. The differences between the experimental data sets at lower values of a for the two post 
mounts are unsettling in that the top mount shows pitch instability near an angle of attack of 0.° The bottom 
mount experimental data and the computations agree better. 

More recently, CFD results have been obtained with the inclusion of the large-post mounting system of the 
14- x 22-ft ST and with the blade-mounted system of the NTF. Starting with the mounting post effects in the 14- x 
22-ft ST, figure 9 shows a plot of C m vs. a computed using FUN3D, albeit with a different, coarser grid than used 
in figure 7. Computed results are shown with and without the post, but in both cases the effect of tunnel walls are 
not included. Clearly the inclusion of the post has shifted the moment coefficient in the right direction, and by the 
right amount. The computed moment coefficient (both with and without post) shows a much broader plateau than 
the experiment around a =15°, or, for that matter, broader than the previous computation without post on the 4.9 
million node grid; this is believed to be a consequence of the relatively coarse mesh. Simulations with a finer 
mesh are currently underway, as are computations including both the post and the tunnel walls. 
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Similarly, interference effects were examined for the NTF mounting system. In this case, only the blade was 
modeled (the part of the support visible in Figure 4), not the large arc-sector located further downstream that 
supports the blade. Furthermore the tunnel walls were not modeled. Figure 10 shows the results computed using 
FUN3D with and without the blade, along with the NTF data. Clearly the inclusion of the blade moves the CFD 
results into better agreement with the data, though, the stability level prior to the lift break (a = 8° or 9°) still 
appears to be over estimated. 

All the results presented above were for low subsonic Mach numbers. Transonic flows have also been 
considered. Figure 11 shows comparison of CFD results, including the NTF mounting blade, computed using 
FUN3D and the unstructured solver USM3D 16 , at conditions corresponding to an NTF test at Mach 0.85 and a Re c 
= 10 million. Results from both codes were generated using the same mesh; note however that since USM3D is a 
cell-based code, the effective resolution of the mesh is approximately six times greater for USM3D than for 
FUN3D. Both codes capture the stability level within the linear range of pitching moment very well, while 
FUN3D predicts an early (relative to experiment) pitch break near a ~ 5°, and USM3D predicts a late pitch break 
at a ~ 6°. 

A final example of static computations are shown in figure 12, which is an example of computations for the 
BWB model tested in the 14- x 22-ft ST in sideslip. For this case, the mounting support was not included in the 
CFD computation. As seen in the results, the predicted computational slope of rolling moment coefficient versus 
sideslip is in good agreement with experiment for values of (3 > 5°. Additional analysis would be required to 
determine if sting effects have influenced the rolling moments for (3 ~ 0°. 


2. Dynamic computations 

Pitch damping for the BWB was chosen as the highest priority dynamic derivative to be computed because of 
the criticality of longitudinal motions for flying-wing configurations. Unfortunately, the national capability to 
measure pitch damping in wind tunnels is very limited, and is usually conducted at low Reynolds and Mach 
numbers. Consequently, the team embarked on time-accurate computations simulating forced oscillations in the 
longitudinal plane, where the model is oscillated at a fixed sinusoidal frequency about a mean angle of attack. The 
resulting variation of pitching moment in this condition quantifies the aerodynamic pitch damping derivative, 
which is a combination of terms due to pitch rate and rate of change of angle of attack. 

Load limitations on the wind-tunnel forced oscillation rig used to obtain the experimental data meant that the 
tests had to be conducted at a low dynamic pressure corresponding to Re c = 700,000. Data were taken in the 14- x 
22-ft ST over a range of mean angles and for several values of pitch amplitude and reduced frequency. For the 
range of mean angles considered here, the corresponding experimental data were taken with the model mounted 
from the bottom. The majority of the computations were carried out for pitch amplitudes of plus or minus 5°, and 
for a reduced frequency, k, of 0.07. None of the force-oscillation computations attempted to date have included 
the post or tunnel walls. For these data curves, the S&C analyst is interested in the overall character of the data 
trend. If the trend of pitching moment with a oscillations results in a counter-clockwise variation of pitching 
moment, the configuration is damped. If, however, the loop exhibits a clockwise trend, the configuration will be 
undamped (propelling) for that condition. The analyst is particularly interested in the magnitude of pitching 
moment at the time of positive and negative pitch motions at the mean value of a since this defines the magnitude 
of the pitch damping derivative for that particular value of a. It should be noted that for all of the plots to be 
shown in this section the loops are counter-clockwise yielding pitch damping with the exception of one plot, 
which will be discussed further. 

As these time-accurate computations were attempted, it became apparent that a real challenge was to 
determine the appropriate balance between grid resolution and temporal resolution for this class of problems. As 
the grid resolution becomes finer, the time step requirement for temporal convergence can become more 
demanding as well. Clearly, even with the increasing availability of computer resources, there is a prudent 
engineering limit to the amount of resources to devote. None of the results obtained to date are believed to be grid 
converged in a strict sense although they are believed adequate from an engineering sense. The time step 
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refinement studies suggest that the results presented below are reasonably well converged temporally for the grids 
used. As for the static moment cases, several codes were employed: COBALT, FUN3D, and USM3D. Further 
details of the FUN3D computations may be found in reference 17, while reference 18 presents the COBALT 
computations in greater detail. In particular both references discuss mesh and time step studies that were used to 
determine practical mesh and time step ranges for subsequent computations with each code. 

Figure 13 shows a comparison of the simulated forced-oscillation results computed with the three codes for a 
mean angle of 0°. USM3D and FUN3D used the Spalart-Allmaras (SA) model while COBALT utilized Menter’s 
(SST) model. 19 The shapes of all three computed curves are similar to the wind-tunnel data in some respects and 
different in other respects; it is difficult to judge which is “better” overall. 

Both the COBALT and FUN3D codes have been used to obtain solutions for mean values of a of 0°, 8°, 16°, 
24°, and 32°. Figure 14 shows a summary plot for FUN3D wherein the computed dynamic data for various mean 
angles are shown in different colors on top of the computed static pitching moment curve, with a similar treatment 
of the experimental data. For the pitch amplitude and reduced frequency considered, the general shape of the 
dynamic pitching moment curve is heavily influenced by the local shape of the static moment curve— which is, in 
turn, shaped by the flow separation states. Although not shown, the corresponding plot for the COBALT results 
would show a very similar trend. 

Detailed results from both codes, for mean angles of 8° and 24° are shown in figures 15 and 16 respectively. 
Generally speaking, the shape of the computed dynamic pitching moment curves agree well with the data for a 
mean angle of 8°, but the agreement has degraded considerably by 16°, as seen in figure 14 for the loop centered at 
a = 16°. Note that between a = 14° and 22°, the experimental static curve in figure 14 exhibits an essentially 
constant static pitching moment, while the computed static pitching moment (from both codes) has a slight 
negative slope. At a mean angle of 24°, any semblance of the computed data to the experimental data has 
vanished, as seen in figure 16. The experimental curve has collapsed near the mean angle, indicating a degradation 
of pitch damping at that value of a. Both computed RANS results show wide loops without significant loss of 
pitch damping. Note also that the experimental static pitching moment coefficient in figure 14 over the range of 
22° to 28° has a much steeper negative slope than the computed static curve. At a = 32°, however, the decrease in 
pitch damping has disappeared from the experimental data and reasonable agreement between computation and 
experiment is once again observed. It is also observed that beyond 30°, the computed and experimental static C m 
curves have nearly the same negative slope. 

To further investigate the poor comparison between experiment and data at 24°, COBALT was used to 
perform computations using DES. 20 Extensive experience with the DES model in COBALT suggested that in 
many high-angle of attack situations, DES, which has inherently less dissipation in off-body separated regions 
than RANS, leads to better agreement with measured data. Since for DES (in the LES region) modeling errors 
decrease as the grid is refined, the original mesh was adapted (based on vorticity) to cluster points in the off-body 
separation region. The SA-based DES results on the adapted mesh are shown in figure 17, as are the RANS results 
on the original mesh (i.e. the same results shown in figure 16), and DES results on the original mesh. In the 
interest of clarity, phase average data of the DES computations are shown, rather than the complete instantaneous 
data set. It should be noted that the DES calculations, while on a per-time-step basis are no more expensive to 
compute than RANS on the same mesh, require perhaps an order of magnitude more time steps before meaningful 
phase-averaged data can be obtained. The fact that the phase-averaged DES data shown figure 17 still have some 
high frequency content indicate that more times steps are really needed before computing the phase averages; 
nonetheless, the trend that DES is giving results more representative of the experiment is unmistakable. Also in 
the interest of clarity, SST RANS results on the adapted mesh are not shown. Only the DES results on the adapted 
mesh show a decrease in pitch damping; the RANS computations on the same adapted mesh do not. 

As an additional note to the discussion about the DES results given in figure 17, the experimental wind-tunnel 
data also have frequency content that has to be smoothed by averaging 40 or 50 cycles of data. Consequently, the 
high-frequency content of the DES solution may be physical and realistic. If the comparisons of figure 17 are 
representative, expending the resources to do the loop averaging for DES may be a requirement, even if it is an 
expensive requirement at this time. 
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3. Summary 


While the computational effort at Langley has been underway for a relatively short time, a good deal of 
progress has been made in advancing the understanding of the requirements for using CFD as a tool for S&C. 
When the work began, the dramatic influence of the support interference on the close-coupled BWB configuration 
was not fully appreciated. However, the results obtained to date suggest that if the support is modeled by the 
computation, then the agreement of the computed static pitching moment coefficient with the experimental data 
can be improved considerably. Still to be investigated is whether or not modeling the support when performing 
dynamic simulations (considerably more difficult and expensive) will yield better agreement with the measured 
data. 

Determining grid and time-step convergence for unsteady flows is extremely resource intensive. Relatively 
speaking, it is far easier to investigate time-step convergence since a halving of the time step simply doubles the 
compute time and (in the asymptotic region) leads to a quartering of the temporal error for the typical second 
order schemes used. In contrast, to quarter the error in the second-order spatial discretization used in most 
’’production" CFD codes requires that the mesh size be halved in each direction, which leads, in the case of 
structured grids, to a factor of 8 increase in the computing time. Compounding the problem, finer meshes 
typically require smaller time steps for temporal convergence, so one is easily faced with a factor of 16 increase in 
resource requirements to reduce the error by a factor of four. In this sense, grid convergence is the tall pole in tent, 
and as mentioned, it is unlikely that the results for the BWB obtained to date are grid converged. Nevertheless, it 
is always important to address grid convergence and it will be discussed further in the status report from AFRL. 
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B. Naval Air Systems Command (NAVAIR) 

The nature of today’s Defense systems acquisition process demands the procurement of increasingly complex 
systems in shorter time periods with less overall funding. In particular, the test and evaluation (T&E) community 
is increasingly challenged to reduce overall test time yet not diminish their mission of comprehensive independent 
assessment of the warfighter’s weaponry. This challenge becomes increasingly difficult without the ability to 
employ advanced technologies aimed at assisting the tester in performing their mission. 

The most promising approach to support this acquisition challenge is virtual testing through the use of 
Modeling and Simulation (M&S). The benefits of M&S are not new, but model validation, computing access and 
speed, and efficient tester interface with multi-disciplinary M&S systems are enabling M&S to more effectively 
support testing programs. 

Under the auspices of the DoD High Performance Computing (HPC) Modernization Office, Naval Air 
Systems Command (NAVAIR) is managing a tri-service software development program that is aimed at 
integrating multiple disciplines through the linking of various codes under a single framework called STEPNET 
(Simulation, Test and Evaluation Process NETwork). These codes and respective application models are then to 
be employed by the engineer to solve any one of a set of aeromechanical problems. The umbrella for the current 
work is the Collaborative Simulation and Testing (CST) Portfolio. CST is an integrated portfolio of high 
performance computing codes directly and specifically aimed at bridging the gap between high-end M&S 
capabilities and test program needs. It is focused on “applied” high-fidelity M&S that makes advancements in 
model validation, creates the ability for data to be quickly processed for “near-time” (as opposed to “real-time”, or 
“near real-time”) support of test programs, and enables effective communication interfaces between test engineers 
and the high-fidelity M&S front-end. Funding for CST comes from the DoD HPC Modernization Program’s 
Common High Performance Software Support Initiative (CHSSI) and is managed by Dr. Dave Findlay. 
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Computational Fluid Dynamics (CFD), wind tunnel tests, flight dynamics analysis tools, and flight simulation 
programs have been used separately and independently of the T&E process in the past. The goal of the sub-project 
"Integrated Simulation of Air Vehicle Stability and Control for Test and Evaluation" is to integrate CFD analysis 
with flight simulation to demonstrate an improved T&E process for fixed-wing aircraft. CFD methodology, based 
on the NASA Langley-developed USM3D CFD code (part of an existing state-of-the-art computational system, 
known as TetrUSS 21 ) in conjunction with NAVAIR-developed CASTLE flight simulation code, will be used for 
the flight test planning and data validation/evaluation process. This project is expected to provide a diagnostic tool 
for flight test decision makers and contribute to the reduction of risks, integration between modeling and 
simulation and testing, and closer coupling of different technical areas. 

The joint NAVAIR/NASA commitment provides a substantial synergy for improving the capabilities of 
TetrUSS, to enhance its ability to solve important problems for stability and control. A primary thrust of this 
project is the efficient computation of dynamic stability derivatives across an aircraft flight envelope for input into 
flight simulation programs. One fundamental requirement is for TetrUSS to be a validated, useable, robust system 
for computing unsteady, often massively separated flows on complex military and civilian aircraft configurations 
in dynamic motion with a known range of uncertainty. Work has been performed to enhance, customize, and 
validate TetrUSS capabilities toward ready computation of dynamic derivatives across the flight envelope. The 
primary support includes the implementation and validation of unsteady Reynolds-Averaged Navier-Stokes 
(RANS) and Detached Eddy Simulation (DES), a forced oscillation and 6-DOF free-to-roll (FTR) analysis 
capability, low-Mach number preconditioning, utilization of advanced two-equation turbulence models, Chimera 
overset moving grids, and assessments of code scalability across multiple HPCMO computer platforms. Work is 
well underway for each of the above capabilities. 22,23 

The examples to be illustrated for this status report are from analyses of the S&C characteristics of the F/A- 
18E/F. This aircraft, in its pre-production version, experienced abrupt wing stall which resulted in its 
susceptibility to asymmetric wing stall, or wing drop. 24 Because of the abrupt wing stall of pre-production 
aircraft, this version makes a good test bed for code comparisons and will be used for the current investigation. 
The wing drop on the pre-production aircraft was eliminated by modifying the leading-edge flap schedule and 
adding a porous surface over the wing-fold fairing. The production F/A-18E/F, which is being delivered to the 
fleet today, does not have wing drop. 

The goal of the present study was to assess the ability of CFD to predict the longitudinal and lateral/directional 
S&C characteristics of the pre-production F/A-18E at transonic Mach numbers without prior knowledge of 
existing wind tunnel or flight test data. Executing the computations in the "blind" would address one of the main 
concerns of the traditional S&C community that the CFD engineer often knew the answer before doing the 
computation while simultaneously educating CFD users about the needs of the S&C community. The study 
summarized in this paper is described fully in a reference by Green. 25 

The geometry evaluated during this computational study was that of the 0.08-scale pre-production F/A-18E 
wind tunnel model. A picture of the computational geometry is shown in figure 18. The wing was modeled with 
67874° flaps, indicating that the leading-edge flap was deflected 6°, the trailing-edge flap was deflected 8° and 
the aileron was deflected 4° downward. The geometry was modeled with both horizontal and vertical tails 
present. Horizontal tail deflections, 6 S , of -6°, 0° and 6° were evaluated for the longitudinal study, while 6 S = 0° 
was used for the lateral/directional study. A Sidewinder missile and launcher were modeled at the wing tip and 
the inlet was modeled as flow-thru. The aft-body distortion required for sting mounting of the 0.08-scale wind 
tunnel model was included in the computational geometry as was the presence of the sting itself. The sting 
geometry was truncated at 50% fuselage length behind the nozzle. 

7. Grid Assessment 

Grids with four different fidelities were used during this study. Coarse, medium and fine grids, see figure 19, 
were used for the predictions of the longitudinal stability and control characteristics. In these three grids, only 
half of the aircraft was modeled and a symmetry-plane boundary condition was implemented. All of the grids 
generated during this study had viscous tetrahedral cells near the surface and inviscid tetrahedral cells away from 
the surface. The y-plus value of the first cell above the surface was approximately unity for each of the grids. 
The coarse grid used 6.5 million cells but had minimal resolution of the verticals tails and the horizontal tail. The 
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medium grid with 11.5 million cells was constructed to bring the resolution of the verticals and horizontal up to 
the level of the wing used in the coarse grid. The fine grid, at 15.1 million cells, started with the medium grid and 
enriched the gridding over the wing. A fourth, full grid for lateral/directional calculations, not shown in figure 19, 
was originally intended to be twice that of the medium grids, or 23 million cells. This, however, turned out to be 
too resource intensive and a compromise between the coarse and medium grids was achieved which yielded a 
total of 16.2 million cells. The lateral/directional grid was generated with 6 S = 0°. 

Since the CFD calculations were to be performed without any prior knowledge of existing wind tunnel or 
flight test data, a grid resolution study was performed to determine the size of the grid necessary to accurately 
calculate the forces and moments on the pre-production F/A-18E. Accurately predicting the longitudinal 
characteristics of the aircraft was the main objective of performing the grid refinement study. In this section, the 
results of the grid refinement study at = 0.8 and 0.9 are presented in figure 20, where lift and pitching moment 
coefficients are presented as a function of angle of attack for the coarse, medium, and fine grids. 

For the = 0.8 results, the coarse grid results are very similar to those of the medium and fine grid results, 
although the pitching moment coefficient of the coarse grid near a = 10° is slightly more negative than that of the 
medium and fine grids. Since values of 6 S = -6° and 6° were going to be evaluated, it was anticipated that this 
difference in pitching moment could become larger for the non-zero tail deflections. As a result, the medium grid 
was used for the calculations with 6 S = -6° and 6° at = 0.8. The sudden slope change in the lift and pitching 
moment curves in figure 20 at a = 6° are indicative of the abrupt stall phenomenon that the pre-production F/A- 
1 8E experienced. 

For the = 0.9 results, see figure 20 again, the lift and pitching moment coefficients for the coarse, medium 
and fine grids are very similar. However, the onset of stall for the fine grid is 1 ° lower in a than that for the 
coarse and medium grids. As a result, there is a significant difference in the lift and pitching moment coefficients 
at a = 7° resulting in earlier instability (positive slope of pitching moment vs. a curve). Without prior knowledge 
of the wind tunnel data, it was assumed that the solutions generated with more points on the wing would be more 
accurate than those solutions with less resolution on the wing (even though wind tunnel data would later suggest 
otherwise). As a result, the fine grid was used to analyze the configurations with 6 S = -6° and 6° at = 0.9. It is 
understood that even with the fine grid, the solutions cannot be considered grid resolved. 

2. Longitudinal Stability and Control 

The longitudinal stability and control predictions for = 0.8 are compared to wind tunnel data from the 
NASA Langley Research Center in Hampton, VA, NASA Ames Research Center in Moffett Field, CA and 
CALSPAN in Buffalo, NY. (More details are again given in the paper by Green. 25 ) In figure 21 both lift 
coefficient and pitching moment coefficients are summarized for the 3 horizontal tails settings. While all three 
experimental entries generated data for 6 S = 0,° only the Ames entry generated data for the non-zero deflection 
angles. Consequently, the experimental data for 6 S = 0° gives insight into the differences between wind tunnel 
entries in different tunnels and are an indicator of the uncertainty in the wind tunnel data due facility and 
installation differences. The CFD results compare well with the wind tunnel data between -2° < a < 6°. 
However, the CFD appears to predict an earlier stall onset at a = 7° than the wind tunnel data indicates. Between 
8° < a < 12°, CFD predicts an approximately linear variation of lift whereas the wind tunnel indicates otherwise. 
The reasons for these discrepancies are not known, but may be related to the unsteady nature of the abrupt wing 
phenomena. Importantly from an S&C perspective, the trends and increments for the horizontal tail deflections 
are well represented in the computations. 

The corresponding pitching moment coefficients are also presented in figure 21 and tell a similar story. First, 
stall seems to be reflected at a lower value of a than in the wind tunnel. However, increments and general shapes 
of the curves for the different horizontal tail deflections are well represented. Green 25 does note that the code 
predicts less stable longitudinal stability characteristics between -2° < a < 6°. 
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3. Lateral/Directional Stability Predictions 

The lateral/directional predictions were performed with the appropriate full-span grid of 16.2 million cells. As 
mentioned already, this full-span grid is more coarse than the medium grid used for the longitudinal predictions. 
In figures 22 to 24, the rolling moment, yawing moment and side force coefficients for the CFD calculations and 
the wind tunnel data are plotted as a function of |3 for = 0.9 for values of a = 5°, 7°, and 10°. The values of (3 
were taken to be -5°, -2°, 0.1° and 5°, where the value of (3 = 0.1° was used in the place of (3 = 0° with the intent 
that any tendencies toward significant flow asymmetry might be identified by the computations. 

The results for a = 5° are shown in figure 22 and highlight the good prediction of rolling moment, yawing 
moment, and side force for these conditions. The paper by Green 25 also includes results at a = 5° and 10° for 
= 0.8 and both of those plots show comparable agreement to what is shown here in figure 22. However, the story 
at this higher value of = 0.9 becomes more complicated than the story for = 0.8. Figure 23 illustrates 
some of the difficulty from the wind tunnel side of the story for a value of a = 7°. Here the data are from both a 
NASA facility and an industrial tunnel and it is obvious that there are significant differences between two 
experimental data sets as conditions approach stall at this higher value of M„. Interestingly, the computational 
predictions, again done in the blind, fall between the two experimental curves. 

The results for a = 10° are shown in figure 24. Here, for the first time, was a real break-down in the blind 
predictions, illustrated by the red symbols. Everything had been predicted so well by the computations up to this 
point that the research team was surprised that the data miss was so dramatic. However, Green 25 was able to 
discern the reason for the mismatch by varying the initial conditions of the computation. For the original ’’blind” 
predictions, the initial conditions were just the free-stream conditions corresponding to the desired value of (3. 
The later runs, conducted after seeing the discrepancies, assumed initial conditions corresponding to the 
converged solution for (3 = 0.1 °. 

The computations may have rediscovered what the experimental world has had to live with for many years. 
That is, when the flow can exhibit hysteresis, or multiple states, active steps have had to be taken in the wind 
tunnel to either quantify or minimize the hysteresis. Often, the solution to minimize hysteresis in the tunnel would 
be to start taking sideslip data from a nominal (3 = 0° initial condition and then move in either the positive or 
negative (3 directions. In this way, the flow states are going from less separation to more separation. This 
approach is similar to what Green has done here by using a nominally symmetric solution as the initial condition. 
More research is required here to determine if the computational differences are truly the result of hysteresis or 
not. Another approach for the CFD might be to run in a time-accurate mode and to see if that would solve, or 
capture, the mismatch for a = 10.° The unsteady nature of the flow over the F/A-18E/F at these conditions is well 
documented. 26,27 

4. Summary 

Under the CHSSI program and the CST portfolio, NAVAIR researchers have taken a fighter configuration, 
the pre-production F/A-18E/F, which is known to exhibit an abrupt wing stall nature, and have generally made 
blind predictions of both longitudinal and lateral/directional characteristics. In general, the TETRUSS code used 
has successfully captured longitudinal trends with horizontal tail increments and lateral/directional increments 
with sideslip for the undeflected horizontal tail case. The shortcomings identified during this investigation include 

(1) the apparently premature prediction of the wing stall when compared to the experimental wind tunnel data and 

(2) sensitivity to the initial conditions when computing sideslip values for = 0.9. The premature prediction of 
wing stall may be due, in part, to the unsteady nature of the actual flow while the code is not being run in a time- 
accurate mode. The significant impact of the initial conditions for the computations should not have been 
surprising because of the similar need to develop processes in the tunnels to control the impact of hysteresis. 

C. Boeing Commercial Airplanes 

Computational fluid dynamics has been used very successfully to develop the high-speed lines of Boeing’s 
commercial transports since the 767. For attached flows, CFD has been able to predict the performance of a 
configuration with a great deal of confidence. This application, however, involves a very small part of the 
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operational flight envelope of an airplane as can be seen in figure 25. There is a need to apply these CFD tools 
throughout the flight envelope because, as seen in figure 26, the larger amount of time currently spent during wind 
tunnel testing targets high-lift lines, loads, S&C, and the flight simulator needs. While great reductions in tunnel 
time required for defining high-speed lines have been made, there has been very little corresponding progress in 
reducing the required time for the other needs. In order to achieve the goals of significant reductions in 
development costs and time to market, the use of CFD must be extended to the corners of the envelope, both flaps 
up and down. 

1. Transonic Pitch Characteristics 

Pitch characteristics for a tail-off transonic transport model are presented in figure 27. At first glance, the lift 
coefficient is nearly identical between CFL3D and wind tunnel test data for the range of a shown. However, a 
close inspection reveals that the CFD results do not fully capture non-linear C Lcc effects in the mid-cx region. 
Pitching moment comparisons reveal similar discrepancies between experimental data and CFD results. General 
trends with a are clearly captured, but CFD is not capturing all of the non-linear C Ma trends in the mid-a region. 
Comparisons at higher values of a reveal that CFD captures the general shape with a bit more instability than 
indicated by the tail-off test data. It is also worth noting that the offset between CFD and wind tunnel in the 
pitching moment curves is generally not of prime concern for commercial transport S&C. The offset may be the 
result of experimental transition strips, model construction or residuals from mounting system corrections, and is 
typically corrected for by trimming with the stabilizer on a conventional commercial transport. It should be noted 
that while the agreement shown in figure 27 is much better than much that has been presented in the previous 
figures, there is always room for improvement. 

2. Spoiler Analysis 

Spoiler reversal CFD predictions are presented in figures 28 for a nacelles-off transonic transport. As 
detailed by Wilkinson, 28 spoiler reversal is caused by adverse shock movement at higher values of a. CFD 
predictions show that at low values of a the shock moves forward as the spoiler deploys, reducing lift and rolling 
the airplane as intended. At higher values of a, however, results show the opposite trend, with the shock moving 
aft as the spoiler is deployed with a consequent increase in lift and adverse rolling moment. Experimental 
comparisons reveal the shock position from CFD is correctly predicted in each case. The importance of 
employing CFD for this case is that the CFD results aided in development of a theory explaining the spoiler 
reversal. Flow analysis revealed that at low values of a the flow past the stowed spoiler is subsonic, while at 
higher values of a the flow past the stowed spoiler is supersonic. 

Subsequent spoiler analysis was extended to a typical transonic twinjet. Figure 29 illustrates slices of total 
pressure in the freestream, highlighting the vortices and wakes in the CFD simulation with outboard spoilers 
deployed (spoilers shown in red). The interaction of the wake and vortex system shedding from the wing on to 
the tail is of particular interest. These computations were performed with CFL3D and include a system of vortex 
generators as indicated in figure 30 and detailed by Tinoco. 29 Figure 30 also shows near- surface streamlines in the 
vicinity of the deflected spoiler and a total pressure slice parallel to the leading edge illustrating the vortices shed 
by the vortex generators. Pressure distributions are compared to test data at inboard and outboard locations for one 
of the solutions in figure 3 1 . Comparisons reveal the solutions are relatively accurate with good correlation on the 
wing surface aft of the spoiler. The outboard location reveals a spike in the CFD data caused by the presence of a 
vortex generator. The distribution and density of pressure ports in the experimental data was inadequate to 
resolve this detail highlighted by the computation. 

3. Assessing the Effect of Reynolds Number on Longitudinal Control 

Reynolds number scaling effects on elevator control authority at a low Mach number are presented in figure 
32, and detailed by Tinoco 29 . CFD results show good correlation with experimental trends and a general increase 
in control authority with increasing Reynolds number out to flight conditions. Figure 33 illustrates that the flow 
physics behind the trend is the decreased flow separation aft of the hingeline as indicated by the reduction in red 
streamlines propagating from the trailing edge as Reynold number increases. 

4. Directional Stability 

Transonic directional stability characteristics were analyzed on a typical twinjet in sideslip using an 
unstructured grid and the CFD++ Navier-Stokes code 30 with the k-e-R t turbulence model. 31 Surface pressure 
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distributions are shown in figure 34 for a transonic condition at positive (3. Shock position differences (denoted 
by the red/green boundary) between port and starboard wings detail the shock movement induced by the sideslip. 
On the starboard wing, the shock position is strengthened and moves aft as the sideslip reduces the effective wing 
sweep for this largely attached flow condition. On the port wing the wing sweep is effectively increased, and the 
shock strength is reduced and moves forward. 


Directional stability (C Np ) comparisons with test data are shown in figure 35. The figure supports 2 primary 
conclusions. First, CFD calculations under predict the level of directional stability. Second, the trend is 
predicted relatively well out to = 0.94. Longitudinal CFD calculations with the same grid provide essential 
insight towards understanding the level shift between the two curves. The lift curve slope (C La ) is a little higher 
than experimental data. This causes the wing/body contribution to directional stability to be a little higher as well. 
Since the wing/body contribution to directional stability is negative, the CFD results end up a little lower than 
experimental data. Figure 36 shows conclusively that the discrepancy in C Np is in the wing/body term by showing 
the fin increment matches wind tunnel data quite well. The fin contribution is obtained by calculating (AC Np ) F in = 
(C Nf / m 0n - (C Nf / m 0ff for both CFD and test data. Results show the fin increment is predicted accurately over a 
range of Mach numbers including near sonic conditions. This is particularly interesting since it shows CFD can 
provide accurate prediction well beyond the typical cruise condition for a conventional transport. It is also 
important to note that these were the first Navier-Stokes calculations performed in sideslip for this analysis and 
may not be grid converged or utilize the most appropriate turbulence model. 

5. Dynamic Derivatives 

Dynamic derivatives predicted by CFD calculations for rotational rate terms are also presented for 3 airplanes 
in figure 37. Results show consistent correlations between predictions using handbook methods for each of the 
configurations and the CFD calculations which were performed at a low transonic condition using the panel 
method A502 (PAN AIR). Since static derivatives (C La , C N(3 ) from the panel code are consistently higher than 
handbook data, it is not surprising that dynamic derivatives are also higher. Consistency in the CFD to handbook 
differences allows use of the data even when the absolute levels are not precisely accurate. 


6. Summary 

By utilizing CFD, Boeing has already saved the company many dollars by reducing the time needed for wind 
tunnel testing. The company is now aggressively attempting to push the application of high-end CFD tools into 
the S&C challenge areas. The current examples describe the use of CFD to provide a better understanding of 
unexpected results from wind tunnel testing, as was the case for the spoiler reversal. High-end CFD, while too 
expensive to employ to populate entire data bases, can provide a level of detail in describing pressures over wings 
that can not be matched by conventional pressure measurement techniques in the tunnel. CFD has also shown the 
ability to predict Reynolds number trends on the horizontal stabilizer at low Mach numbers and has demonstrated 
the ability to predict the effect of M^ on directional stability. Finally, even employing lower order CFD such as 
panel methods to predict dynamic stability characteristics is shown to provide consistent increments from 
configuration to configuration. 

D. Air Force Research Laboratory (AFRL) 

The Computational Sciences Branch and Aerodynamic Configuration Branch of Air Vehicles Directorate, Air 
Force Research Laboratory (AFRL/VAAC and AFRL/VAAA, respectively), have a small effort underway to 
better understand the fundamental reasons CFD typically does a poor job of predicting S&C. The goal is to 
determine what level of modeling fidelity is required to accurately predict stability quantities, i. e. what level of 
geometric fidelity, grid resolution, and turbulence modeling are necessary to get accurate results. The 
presumption is that the Navier-Stokes equations completely and correctly model this system, and it is just a matter 
of proper application of current (2 nd order) CFD tools that is necessary. The approach is to apply a very 
methodical analysis process to a relatively simple stability problem, addressing each of the areas in turn to 
determine best practices for applying CFD to S&C problems. Therefore, a relatively simple geometry that has 
been tested experimentally, and yet which exhibits somewhat complex behavior in its stability characteristics, was 
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considered. To simplify even further, initially only static pitching moment was considered, as even this most 
basic stability quantity has proven difficult to accurately predict with CFD. 

The 1303 generic UCAV geometry was chosen because there is recent experimental data available at benign 
test conditions, it has a relatively simple delta wing planform, and it exhibits nonlinear behavior in pitching 
moment with angle of attack. 32 There has also been some effort to ensure that the digital representation of the 
model is consistent with what was tested in the wind tunnel. The 1303 model was tested with various leading 
edge geometries, ranging from sharp to varying levels of roundness. Figure 38 shows the model setup during 
wind tunnel testing, while figure 39 indicates the pitch break at moderate angles of attack in the experimental data. 
Surface oil flow patterns, not shown, on the model indicated there is some type of leading edge rollup or vortex as 
the angle of attack increases, which likely results in the observed pitching moment behavior. 

The sharp leading edge geometry was chosen for comparison, the intent being to eliminate possible issues 
with transition over a rounded leading edge. It is apparent from the model setup in the tunnel that the mounting 
sting and apparatus is large relative to the model. For that reason the initial geometry considered with CFD 
includes a simplified representation of the sting assembly (see figure 40). The portion of the sting close to the 
model is an accurate representation of the tunnel apparatus, but after a certain distance the sting cross section is 
held constant and extended to the edge of the computational domain. No other tunnel geometry was considered, 
so each of the following cases was run with freestream conditions at the farfield boundary. A number of 
increasingly refined unstructured meshes were created with VGRIDns, with cell counts of 566K, 1.7M, 4.7M, and 
8.3M cells (see figure 41). The farfield boundary was a box approximately 17 body lengths away from the UCAV 
in each direction; symmetry conditions were assumed to reduce the grid sizes. The smallest grid was an Euler 
tetrahedral grid, while the others were viscous grids with prismatic cells close to the body. Anisotropic stretching 
of 25 to 1 was used for the cells at the leading and trailing edges of the model. The larger grids increased the grid 
density primarily about the leading and trailing edges of the model. The computations were run to match the test 
conditions of a Mach number of 0.25, and Reynolds number of 4.28 million. For all viscous grids this resulted in 
an average y + of 0.5. All computations were run with the Air Vehicles Unstructured Solver (AVUS), formerly 
known as Cobalt 6 o H . The viscous cases were run fully turbulent with either the Spalart-Allmaras one-equation 
turbulence 33 model or Menter’s SST two-equation turbulence model. 34 

Figure 42 compares the computed drag polar to the corrected wind tunnel data. As expected, the inviscid 
solution is the farthest off the data, but it still does a reasonable job predicting the correct character of the solution. 
At the lower values of drag both turbulence models do a reasonable job on both grids. The turbulent results on the 
1.7M cell grid depart noticeably from the experimental data at the higher drag values, while the 4.7M cell grid 
tracks the data well. The choice of turbulence model has little impact on the results for this larger grid. So, as 
expected, CFD predictions of lift and drag when the flow is predominantly attached are generally quite good. 

Figure 43 compares the computed pitching moment coefficient to both the corrected and uncorrected wind 
tunnel data. Since the computations do not include the tunnel walls, the expectation is that the corrected tunnel 
data would be best for comparison. As expected the inviscid results are quite far off, but still reproduce the 
character of the data to a relatively high angle of attack. The viscous results on the 1.7M cell grid mimic the 
inviscid results, indicating the grid does not have the necessary resolution to accurately predict the viscous 
character of the flow (regardless of turbulence model). The 4.7M cell grid corrects this, and for lower angles of 
attack tracks the corrected tunnel data fairly well. The turbulence models show slight differences. At the pitch 
break the computed results depart somewhat from the data, and maintain a somewhat constant offset from the data 
through the last computed angle of attack. Finally, one data point was computed on the largest 8.3M cell grid 
slightly after the pitch break, and shows a noticeable improvement relative to the 4.7M cell grid. Thus, the 
indication is that a sufficiently resolved grid would match the experimental data. But, an 8+M cell grid for such a 
simple planform is much larger than experience with lift and drag coefficients would suggest. 

Figure 44 shows computed surface oil flow patterns and pressure contours for two values of angle of attack 
(4° and 8°) on the 4.7M cell grid with Menter’s SST turbulence model. The solution at the higher angle of attack 
clearly indicates an attachment line for a leading edge roll-up/vortex, indicating the solution has captured the 
character of the experimental results. Figures 45 and 46 show comparisons of computed pressure coefficient with 
experiment at the 18-inch span location for these two angles of attack (semi-span is 30 inches). Surprisingly, at 
the lower angle of attack a major difference between the turbulence models is apparent, with Spalart-Allmaras 
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showing much better agreement with the data on the least resolved grid. One would assume that this agreement is 
fortuitous. The higher angle of attack shows better agreement between both turbulence models and the data, 
though the two models under- or over-predict the data in different locations. In any case, the pressure coefficient 
comparisons show that even though the integrated quantity agreement improved with higher grid resolution, the 
discreet data agreement does not necessarily improve everywhere. 

The next step in the process is to apply solution adaptation techniques to refine the grid only in the regions that 
need it. This appears to be the only way to obtain the accuracy that is needed, while (hopefully) keeping the grid 
size manageable. It may be necessary to add further geometric fidelity, such as the tunnel walls and sting 
assembly. As the flow separates at higher angles of attack, turbulence model differences will likely become more 
pronounced. This may then require that more advanced turbulence modeling techniques, such as detached eddy 
simulation (DES), be applied. 

7. Summary 

The effort is aimed at examining the grid requirements and turbulence model needs to achieve accurate results. 
This UCAV configuration has been chosen for this study because of its simplicity. It has been demonstrated that 
achieving convergence with more grid resolution in an integrated sense of forces and moments do not insure, as is 
well known, that individual pressure distributions are also converging. 
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Summary 

This paper documents the status of the efforts at four organizations to aggressively apply high-end CFD codes 
to the real-world problems and challenges associated with the prediction of the stability and control (S&C) 
characteristics of a new flight vehicle. The informal working relationship between NASA Fangley, Naval Air 
Systems Command (NAVAIR), Boeing Commercial Airplanes (BCA), and the Air Force Research Faboratory 
(AFRF) has been helpful in general and, particularly, because it allows all four organizations to have insights into 
the issues involving the range of vehicles being studied by these groups. These vehicles include, respectively, the 
Blended Wing Body (BWB), the F/A-18E Super Hornet, commercial transport configurations, and a notional 
Uninhabited Combat Air Vehicle, or UCAV. Fessons learned can be summarized for each organization. 

The work by NASA Fangley involved both static characteristics as well as ground-breaking efforts in the 
evaluation of BWB pitch damping. While the work is still underway, it is clear that further improvements are 
needed before consistent prediction of static longitudinal stability for the low-a region are possible. Impacts of 
wind-tunnel blade and post supports are predicted by the codes and appear to do a good job of estimating the 
offsets. In the area of pitch damping, it is found that RANS codes do a credible job of predicting pitch damping 
trends for lower values of a but may be missing the reduction in pitch damping measured in the tunnel near a = 
24°. DES, which may be better equipped to track the complex separated flows through less dissipation, may be 
showing some advantage for that condition. Again, further work is needed in this area. 

One of the objectives of the work at NAVAIR with the F/A-18E has been to address the concern of the 
traditional, experimentally driven S&C community that many of the computational predictions have been made in 
the past with the CFD analysts already knowing the answer. Consequently, the NAVAIR work was carefully 
executed to make sure that the researcher did not know the answer before conducting the grid convergence studies 
or choosing any other running options. The ’’blind” longitudinal results represented trends in lift and pitching 
moment satisfactorily but disagreed with wind tunnel data for angle of attack at which wing stall occurred. The 
lateral ’’blind” results shown were particularly good for predicting rolling moment as a function of sideslip for a 
Mach number of 0.8. The results at a Mach number of 0.9 were mixed with good agreement at lower values of 
angle of attack but completely missing the mark for an angle of attack of 10.° Computations done after the initial 
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mismatch, and not ’’blind,” demonstrated that the initial conditions for the flow solution are important and give a 
very representative answer if the initial condition is taken as a computational solution appropriate for no sideslip. 
In retrospect, this is similar to the wind tunnel, where processes have been developed to minimize hysteresis 
effects in sideslip. 

The work at BCA represents the most applied usage of the high-end CFD because that organization has been 
engaged on applying the CFD to that company’s real-world challenges to augment and complement information 
that they develop in the wind tunnel. From the examples included in this paper, it is clear that the use of CFD was 
extremely useful for the spoiler reversal challenge. The codes not only confirmed what was being seen in thewind 
tunnel but also gave the critical insight to explain the flow physics. High-end CFD has also been used to 
successfully capture both Reynolds number trends and Mach number trends. Consistency in correlation across 
multiple airplanes is critical. Even panel methods can be employed incrementally when discrepancies are 
consistent across multiple airplanes as shown in the linear range dynamic derivative example. 

The effort at AFRL is focused on learning about grid resolution and turbulence modeling requirements 
necessary to achieve an accurate solution. This is work still in progress but has already reiterated that it is 
necessary to examine pressures as well as integrated forces and moments to insure that solutions are truly 
converging as grid refinement is taking place. 

This paper has summarized the work within the four organizations represented. The work has not yet achieved 
its goals of validation of the computational tools for the S&C challenges but is continuing so that these valuable, 
complementary tools to the existing wind tunnel tools can be more routinely integrated into the aircraft 
development cycle. With the ability of CFD to simulate (1) flight Reynolds numbers, (2) flight Mach numbers, 
and (3) free-air conditions without tunnel supports or walls and to provide unprecedented flow diagnostics, the 
practical implementation of CFD in the S&C arena will lead to reduced development times and more optimal 
designs. 
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Figures 
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Figure 1. Examples of high performance military aircraft that have experienced transonic, uncommanded 
lateral motions, such as wing rock or drop. 
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Figure 2. Examples of other aircraft having unexpected stability and control characteristics in flight. 
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Figure 3. Legacy differences between the CFD community and the S&C community. 
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Figure 4. Various BWB model installations used during the experimental program. 
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* NTF, Re=3.7M 

a 14x22, Large Post, Re=3.7M 

e 14x22, Large Post, Re=2.6M 

14x22, Small Post. Re=2.6M 

_l I l_ _l _l 
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a, deg. 

Figure 5. Comparison of experimental data sets for static BWB pitching moment; 
Moo = 0.2 (nominal). 



Figure 6. Key regions of interest for S&C. Region I is longitudinally stable and generally represents 
attached flow. Region II contains the unstable pitch break and continuing breakdown of flow. Region III 
represents a stable flow after stall. 
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Figure 7. Comparison of computed static BWB pitching moment with experiment; 
Re c = 3.7 million, = 0.2. 



Figure 8. Comparison of computed static BWB pitching moment with experiment; 
Re c = 0.7 million, = 0.1. 
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Figure 9. Effect of 14- x22-ft ST large post mount on computed static pitching moment for Re c = 3.7 
million, = 0.2. 



Figure 10. Effect of NTF blade mount on computed static pitching moment for 
Re c = 3.7 million, = 0.2. 
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Figure 11. Comparison of computations including the blade support and experimental static pitching 
moment; Re c = 10 million, = 0.85. 



Figure 12. Comparison of computed and experimental rolling moment versus sideslip; a = 0°, Re c = 3.75 
million, = 0.20. 
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a, deg. 


Figure 13. Dynamic pitching moment as a function of a for a forced oscillation of 
+/- 5 degrees about a = 0°. Re c = 0.7 million. 


E 

o 



-10 0 10 20 30 40 

a, deg. 


FUN3D, 1 BOO steps/cycle, 4.9M nodes, SA 

■a 14x22, Bottom Mount 


Figure 14. Dynamic and static pitching moments as a function of a; forced oscillations of +/- 5 degrees; 
Re c = 0.7 million. 
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Figure 15. Dynamic pitching moment as a function of a for a forced oscillation of 
+/- 5 degrees about a = 8°. Re c = 0.7 million. 



Figure 16. Dynamic pitching moment as a function of a for a forced oscillation of 
+/- 5 degrees about a = 24°. Re c = 0.7 million. 
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Figure 17. DES and RANS simulations from COBALT for oscillations of 
+/- 5 degrees about a = 24°. Re c = 0.7 million. 
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Figure 18. Computational geometry of the 0.08-scale wind tunnel model of the pre-production F/A-18E 
Super Hornet. 
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Figure 19. Grids based on F/A-18E wind tunnel model geometry with 6°/8°/4 0 flaps. 
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Figure 20. Grid refinement study. 
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Figure 21. Lift and pitching moment for 3 horizontal tail deflections. = 0.8, Re c = 3.9 million. 
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Figure 22. F/A-18E in sideslip, = 0.9, Re c = 4.0 million, a = 5°. 
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Figure 23. F/A-18E in sideslip, = 0.9, Re c = 4.0 million, a = 7°. 
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Figure 24. F/A-18E in sideslip, = 0.9, Re c = 4.0 million, a = 10°. 
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CFD mostly done 



Figure 25. CFD and the flight envelope. 
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Figure 26. The opportunity to expand CFD usage and reduce wind tunnel testing hours. 
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Figure 27. Pitch characteristics of a tail-off transonic transport. 
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Figure 28. Transonic spoiler reversal. 
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Figure 29. Transonic Navier-Stokes simulation with spoilers, which are in red. 



Figure 30. Flow detail near spoiler of figure 28. 
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Figure 31. Pressure distribution with deflected spoiler. 
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Tail Reynolds Number 


Figure 32. Elevator control effectiveness. 
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Figure 33. Reynolds number effects with deflected elevator. 



Figure 34. Transonic twinjet in sideslip. 
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Figure 35. Directional stability for a transonic twinjet. 



Figure 36. Fin contribution to directional stability. 
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Figure 37. Dynamic derivative correlation between Panel Method (CFD) and handbook. 



Figure 38. 1303 UCAV in 5m pressure tunnel at Farnborough, England, July 2002. 
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Figure 39. Nonlinear pitching moment for 1303 UCAV with various leading edge geometries. 



Figure 40. Initial 1303 geometry for CFD showing representative sting connected to model. 
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Figure 41. Surface triangulations for 1.7, 4.7, and 8.3 million cell grids (respectively) used in mesh 
refinement study. 



Figure 42. AVUS drag polar comparison with corrected experimental data for various grid resolutions and 
turbulence models (M^ = 0.25, Re = 4.28 million). 
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Figure 43. AVUS pitching moment comparison with corrected and uncorrected experimental data for 
various grids resolutions and turbulence models (M x = 0.25, Re = 4.28 million). 



Figure 44. Surface oil flow and pressure contours on 4.7 million cell grid for solutions with Menter’s SST 
turbulence model at 4° and 8° angle of attack, respectively (M^ = 0.25, Re = 4.28 million). 
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Figure 45. Computed pressure coefficient comparison with experiment at 18-inch span location for 4° 
angle of attack (Moo = 0.25, Re = 4.28 million). 
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Figure 46. Computed pressure coefficient comparison with experiment at 18-inch span location for 8° 
angle of attack (Moo = 0.25, Re = 4.28 million). 
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